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Abstract 

The RI/MOM non-perturbative renormalization scheme is studied on 
the lattice in SU(3) quenched QCD with Wilson fermions. The gauge 
dependence of some fermion bilinear renormalization constants is discussed 
by comparing data which have been gauge-fixed in two different realizations 
of the Landau gauge and in a generic covariant gauge. The very good 
agreement between the various sets of results and the theory indicates that 
the numerical uncertainty induced by the lattice gauge-fixing procedure is 
moderate and below the statistical errors. 
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1 Introduction 



Non-perturbative renormalization techniques [|], [§ have become a crucial ingredi- 
ent in lattice computations of fundamental QCD parameters and hadronic matrix 
elements. They allow the limitations of lattice perturbation theory to be over- 
come in the determinations of the relevant renormalization constants by removing 
the uncertainties due to unknown higher order terms. 

The RI/MOM technique proposed in Refs. |], || has been successfully applied 
to compute renormalization constants of composite fermion operators in many 
lattice regularizations: Wilson fermions [l|], P|-[lQ|, Kogut-Susskind fermions 
[ |TT|j , domain- wall fermions and overlap fermions ||I5fl . The method imposes 
RI/MOM renormalization conditions on conveniently defined amputated Green 
functions computed non-perturbatively between off-shell quark states at large 
virtuality in a fixed gauge. In practical computations the renormalization scale 
/i 2 , determined from the virtuality p 2 of the external states, must satisfy the 
condition Aqqtj <C \i <C 0(1 /a) to ensure the reliability of the results [|TJ. 

The RI/MOM renormalization constants depend in general on the external 
states and therefore on the gauge in which the RI/MOM renormalization condi- 
tions have been imposed. For gauge-invariant operators the coefficients needed 
to match in a given gauge-invariant scheme will cancel these dependences up to 
higher orders in the continuum perturbative expansion and up to discretization 
errors. A necessary assumption in this procedure is that, at large virtualities, 
the continuum limit of the renormalized gauge-dependent matrix elements, com- 
puted by imposing the non-perturbative Landau lattice gauge fixing, lead to the 



corresponding values computed in perturbation theory up to higher order [14 



The existence of both continuum [|15[] and lattice Gribov copies fl6| , [lTj and the 
numerical noise that they can generate [|18[] (for a recent review see Ref . |19| ) is an 
unsolved problem of the lattice non-perturbative gauge fixing. The real concern 
is, of course, the influence that these phenomena may have on the values of the 
renormalization constants, when computed using non-gauge-invariant quantities 
as in the RI/MOM scheme. 

The aim of this paper is to study the systematics induced by the gauge-fixing 
procedure on the RI/MOM determinations of the renormalization constants. We 
have computed the quark field (Z^), the axial- vector (Za) and scalar (Z$) renor- 
malization constants, by imposing the standard lattice Landau gauge and another 
realization of the Landau condition that, in the following, will be called Landaul 
gauge. Whereas the two realizations impose the same gauge-fixing condition in 
the naive continuum limit, they are affected by different Gribov ambiguities. By 
comparing the different sets of results for the renormalization constants, we have 
found negligible differences within our statistical errors. 

As will be discussed in section 2, a very interesting feature of the Landaul 
gauge is that it can be generalized to impose a generic covariant gauge on the 
lattice, as proposed in Ref. [O. By exploiting this opportunity, we have extended 
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the study of the gauge dependence to off-shell Green functions measured in a 
generic covariant gauge. 



2 Gauge-Fixing Conditions 



In this section we briefly describe the continuum analogous of the non-perturbative 
gauge-fixing conditions whose discretized forms have been used to obtain the nu- 
merical results presented below. 

In the standard Landau gauge the expectation value of a gauge-dependent 
operator is given by 



(O f ) = J 5A^Oe- s ^/\ F [A]5{d^A 

where the Faddeev-Popov factor is defined as 

A F [A] =det (d^D^A]) 
and the gauge-fixing condition 



d^A G = with periodic boundary conditions . 



is enforced non-perturbatively by minimizing the Gribov functional JLS 

F A [G]^\\A G \\ 2 -- 



d x Tr 



A G A G 



(1) 



(2) 



(3) 



(4) 



Since the functional in Eq. @ can have multiple minima in the non-perturbative 
regime, the corresponding gauge-fixing condition is incomplete, a fact that is 
known as the Gribov copy problem . 

In order to study the effects induced by Gribov's ambiguities on gauge- 
dependent matrix elements, the Landau gauge can also be fixed by minimizing a 



different functional [14 



H A [G] 

which reaches its minima when 



d A xTr 



(MS 



G\2 



(5) 



(6) 



where D u is the covariant derivative. Equation (|6]), supplemented with periodic 
boundary conditions, defines a new gauge which we call Landaul in the following. 
In this case the expectation value of a gauge-dependent operator is given by 



(O h ) = J 5A,Oe- s ^/\ H \A]5(D v d p d,A, 



(7) 



where A#[A] is defined analogously as in Eq. (Q) JHJ. By comparing Eqs. (|1|) 
and (0) it is straightforward to show that the two procedures lead to the same 



values of gauge-dependent matrix elements in the perturbative regime [20]. 

The solutions of Eq. @ are all absolute minima of Ha[G], at variance with the 
copies associated with the standard functional in Eq. (^). Therefore we expect a 
different distribution of Gribov's copies for the two procedures, which could gen- 
erate different effects in off-shell matrix elements. It is clear that Eq. (U) could 
also have "spurious" solutions, which correspond to zero modes of the operator 
D u d u and do not satisfy the gauge condition in Eq. (|]). The numerical results 
of the exploratory study presented in Ref . |L4| , and, a posteriori, the new results 
presented below, indicate that such a possibility can be discarded. 



The functional Ha[G) has been proposed in the literature |2T|, [14]] because it 
can be easily generalized to the following form 



(A,A) 



[G] 



d xTr 



(d»A° - A)(d v A« - A) 



G 



(8) 



where A(x) are matrices belonging to the Lie algebra of the group and are gen- 
erated according to a Gaussian distribution with the variance £ fixed at a given 
value. It is the value of this parameter that determines the corresponding covari- 
ant gauge. The H(a,a) [G] minima satisfy the gauge condition 



Dydy^d^A^ — A) = , with periodic boundary conditions . 
The expectation value of a given operator is obtained by 

(0ir)e = J 5Ke-y dixTl[K2] j6A,Oe- s ^A H [A]5(d,A,~A) 



(9) 



(10) 



For a complete description of the method, see Ref. |L4]]. In the following we 
perform a further check of the gauge-dependence comparing the results obtained 
in the Landau gauge and in one of the lattice covariant gauges. 



3 Numerical Results and Discussion 



We have used a sample of 80 SU(3) gauge configurations retrieved from the repos- 
itory at the "Gauge Connection" flhttp:/ /qcd.nersc.gov/D , which were generated 
with the Wilson gluonic action at /3 = 6.0 and V xT = 16 3 x 32 
discretized version 



sec 



14}] for details) of the functionals in Eqs. 



By using the 
, (|) and (§) 

we have rotated each configuration in the Landau gauge, in the Landaul gauge 
and in the covariant gauge with the bare value of the gauge fixing parameter 
£ = 8 respectively. 

The various gauges have been enforced with a quality factor 9 < 10 



sec 



14 



for details) with the Origin2000 computers at Boston University. The inversion of 
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p 


6.0 


Action 


Wilson 


Lattice Size 


16 3 x 32 


# Configurations 


80 




Landau 

T , £1 T~) H £1 11 1 

£ = 8 


k 


0.1530 
0.1535 
0.1540 
0.1545 
0.1550 
0.1555 


k c ra 


0.15703(2) 


a' 1 [23] 


2.26(5) GeV 



Table 1: Main parameters of our simulations. 



the Wilson-Dirac operator has been performed on an INFN APE100 machine for 
masses corresponding to the values of the k parameter shown in Table [l]. Once 
the quark propagator S(x, 0) = (ip(x)ijj(O)) has been computed for each mass and 
gauge and Fourier-transformed, we have determined the inverse Euclidean bare- 
quark propagator S~ 1 (pa). We have computed the amputated Green functions 
of the local quark bilinears Or{x) = ip(x)rifj(x) in the Fourier space Q 



Ar (pa) = S 1 (pa)Gr(pa)S 1 (pa) 



where 



Gr(p) = / d A xd 4 ye 



-ip(x—y) 



Mz)Or(0My)> 



(11) 
(12) 



and T is a generic Dirac matrix corresponding to T — {A, S} = {7^75, 1}. 

In the RI/MOM scheme, the wave-function renormalization constant is 
defined as 



1 



lim — i — Tr 

ro->0 12 



as- 1 (pa) 



dp 1 



In order to avoid numerical derivatives, it is convenient to use 

1 Tr £,,=1,4 7 M sin(p At a)S'- 1 (pa) 



ZJfia) 



lim — i- 

m->0 12 



4Eu=i, 4 sin (p M a) 



(13) 



(14) 



which differs from Z^ by a finite term of order a s . The matching between the two 
renormalization constants can be computed in continuum perturbation theory. 

2 In what follows we adopt the conventions of Ref. || . 
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Since this correction is negligible at the level of the numerical precision of our 
calculation (see below), we will not take it into account in the following. Therefore 
throughout this paper we assume Z^ — Z'^. 

The renormalization constant Z^ l (p,a, go); which defines the renormalized op- 
erator O^ 1 = Z^O r , is fixed by imposing in the chiral limit the renormalization 
condition 

Z^CHVMTr P r A r (pa)\ p2=ll2 = 1, (15) 
where Pr is a suitable projector on the tree- level amputated Green function 

Even though Oa and Os are gauge-independent operators, their matrix ele- 
ments between quark states acquire a gauge dependence. In Fig. [I] are shown 
the numerical results for Z^, Za and Z$ calculated in the Landau and in the 
Landaul gauge, as a function of the square lattice momenta 

M 2 = 4sin 2 (^). (16) 

In the perturbative region, Z^{a[i) has a flat dependence on the renormalization 
scale. This is expected since in the Landau gauge it has zero anomalous dimension 
at the leading order. Also shows a flat behaviour due to the compatibility 

of the RI/MOM renormalization conditions with the axial Ward identities. The 
value of Zs{a[i) increases with a logarithmic behaviour, as expected in pertur- 
bation theory (see also Fig. [|). Our data are in very good agreement with the 
results reported in the literature 0. 

The results for all correlators, corresponding to Landau and Landaul gauges, 
coincide within the statistical errors. As a consequence we can conclude that the 
data do not show any residual effect due to the presence of lattice Gribov's copies 
in the statistical sample of configurations generated and for the lattice used. This 
represents one of the main results of this paper. 

To show the amount of gauge dependence that can be found in our correlators 
and therefore our sensitiveness to the gauge, in the last part of this section we 
compare the results obtained in the Landau and in the covariant gauge with the 
bare gauge parameter fixed to £ = 8. The gauge dependence of the primary 
quantities Z^, A a and A s measured on the lattice is clearly visible in Fig. 0, 
where they are plotted for the Landau and the covariant gauge as a function 
of the momenta {ap) 2 . All three quantities show a statistically significant gauge 
dependence while the renormalization constants Za, Z$ shown in Fig. |2|, obtained 
by computing the ratios as indicated in Eq. fll5"D, do not exhibit any visible gauge 
dependence within the statistical precision reached in our simulation. For these 
quantities, the fluctuations of the simulation hide the weak gauge dependence 
that is expected in perturbation theory from the next-to-leading order terms in 
a s . These results indicate an upper limit to the numerical troubles that can be 
expected in the RI/MOM renormalization constants because of the lattice gauge 
fixing. 
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Figure 1: The Za and Z$ renormalization constants are shown as a function of 
(an) 2 for the Landau and Landaul gauges, as discussed in the text. The data 
have been slightly displaced in the x-direction and the error bars for one set only 
have been shown in order to help the eye. The error bars for the other sets are 
comparable. The errors are jackknife. 
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Figure 2: The behaviour of Z^, Ka and As is shown as a function of (ap) 2 for 
the Landau (black dots) and the covariant gauge (crosses) data. The errors are 
jackknife. 
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Figure 3: The behaviour of Za and Z$ is shown as a function of (a/z) 2 for the 
Landau (black dots) and the covariant gauge (crosses) data. The errors are 
jackknife. 
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4 Conclusions 



We have proposed a method to study the effects of Gribov's copies on lat- 
tice gauge-dependent matrix elements. We have implemented it numerically in 
quenched QCD with Wilson fermions and we have computed the matrix ele- 
ments needed to fix the quark field, the axial-vector and scalar renormalization 
constants in the RI/MOM non-perturbative renormalization technique. Our nu- 
merical results do not show any significant differences between data coming from 
the standard Landau and the Landaul gauge. We do not find any relevant effect 
that can be traced back to the incompleteness of the gauge fixing for a number 
of gauge-field configurations and for lattice sizes of the order of those used in this 
calculation. The main indication of our paper is therefore that, at the level of the 
precision that we have reached, the systematics due to the numerical gauge-fixing 
machinery has negligible effect on the matrix elements needed for implementing 
the RI/MOM scheme. 
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